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Abstract. With the help of families of density contrast indicators, we study 
the tendency of gravitational systems to become increasingly lumpy with time. 
Depending upon their domain of definition, these indicators could be local or global. 
We make a comparative study of these indicators in the context of inhomogeneous 
cosmological models of Lemaitre— Tolman and Szekeres. In particular, we look at the 
temporal asymptotic behaviour of these indicators and ask under what conditions, 
and for which class of models, they evolve monotonically in time. 

We find that for the case of ever-expanding models, there is a larger class of 
indicators that grow monotonically with time, whereas the corresponding class for 
the recollapsing models is more restricted. Nevertheless, in the absence of decaying 
modes, indicators exist which grow monotonically with time for both ever-expanding 
and recollapsing models simultaneously. On the other hand, no such indicators may 
found which grow monotonically if the decaying modes are allowed to exist. We also 
find the conditions for these indicators to be non-divergent at the initial singularity 
in both models. 

Our results can be of potential relevance for understanding structure formation 
in inhomogeneous settings and in debates regarding gravitational entropy and arrow 
of time. In particular, the spatial dependence of turning points in inhomogeneous 
cosmologies may result in multiple density contrast arrows in recollapsing models over 
certain epochs. We also find that different notions of asymptotic homogenisation may 
be deduced, depending upon the density contrast indicators used. 



1. Introduction 

An important feature of classical self-gravitating systems is that in general they are not 
in equilibrium states. This instability gives rise to spontaneous creation of structure 
which is assumed to lead to increasing lumpiness as time increases. 

This rise of structure with time in self-gravitating systems seems to run counter to 
the usual intuition in classical statistical physics where increasing time is associated 
with the increase in microscopic disorder and hence to macroscopic uniformity in 
matter. The main reason for this difference is thought to be due to the long range 
and unshielded nature of the gravitational force . This dichotomy has led to the 
question of possible existence of gravitational entropy and its connection with the 
usual thermodynamic entropy (see for e.g. |2^, and references therein). 

Whether a satisfactory notion of gravitational entropy exists and whatever its 
nature may be, it is of interest to find out the extent to which the assumption regarding 
the increase in structuration with increasing time holds in relativistic cosmological 
models. This would be of value for a variety of reasons, including its potential relevance 
for debates concerning structure formation in the Universe. Furthermore, the presence 
or absence of indicators that evolve monotonically in time could inform the debates 
regarding gravitational entropy and in particular whether, if it in fact exists, it should 
be monotonic in models which possess recollapsing phases (see for e.g. [[l4| f2l[ p^|). 

Since indicators that best codify such structuration are not known a priori, rather 
than focusing on a single indicator, we shall here consider families of spatially co- 
variant indicators which measure the density contrast and which include as special 
cases indicators put forward by Bonnor ||] and Tavakol & Ellis [^8). We shall also 
consider for completeness some indicators previously used in the literature, including 
that given by Szafron & Wainwright |26| and the non-covariant indicators of Silk 
and Bonnor even though the latter two may, in view of their lack of covariance, 
be viewed as suspect from a physical point of view. 
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We shall employ inhomogeneous cosmological models of Lemaitre-Tolman (LT) 
Ji"9| . p9[ and Szekeres |27| in order to make a comparative study of these indica- 
tors in the general relativistic inhomogeneous settings. Since these models involve 
arbitrary! functions, the integrals involved in the definitions of our indicators cannot 
be performed in general. As a result, we look at the asymptotic behaviour of both 
ever-expanding and recollapsing families of models as well as their behavious near the 
origin, and in particular look for conditions under which the asymptotic evolution of 
these indicators is monotonic with time. Clearly these can only give necessary condi- 
tions for the all-time monotonic evolution of these indicators. To partially extend these 
results to all times we also calculate these indicators for a number of concrete models 
given in the literature. The crucial point is that our asymptotic results bring out some 
general points that seem to be supported by our all-time study of the concrete models. 

The organisation of the paper is as follows. In section 2 we introduce and motivate 
families of density contrast indicators. Sections 3-6 contain the calculations of these 
measures for LT and Szekeres models respectively. In section 7 we consider the be- 
haviour of the dimensionless analogues of these indicators. In section 8 we study the 
behaviour of these indicators near the initial singularity. Section 9 gives a discussion of 
some of the consequences of our results and finally section 10 contains our conclusions. 

Throughout we use units in which c = G = 1 and lower case latin indices take 
values to 3. 



2. Density contrast indicators 

Intuitively one would expect the rise of structuration in cosmology to be related to 
the coarse grained spatial variance of various physical scalars. An important scalar in 
cosmology is the energy density p (or in a multi-component version of this, densities 
pM corresponding to the different components of the content of the Universe) . Here for 
simplicity we confine ourselves to the one component case and consider the rotation- 
free dust setting. We can then introduce a global spatial variability index defined as 



/ 



EzMjv (1) 
Po 

where po is the mean density defined appropriately and dV is the comoving volume 
element. One can make this notion spatially covariant by, for example, expressing it 
in terms of the fractional density gradient introduced by Ellis & Bruni j?J , 

where h a \, = g a b + u a Ub projects orthogonal to the unit 4- velocity u a , which we shall 
throughout assume to be uniquely defined. A related covariant spatial variability 
index can then be defined thus 

\Xa\dV (3) 



f These functions are not in fact totally arbitrary as they need to satisfy certain constraints which 
we shall discuss below. 
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where the integration is over a 3-surface £ or part thereof. 

Now it is not a priori clear what are the indicators that best codify such 
structuration and in particular what their monotonic properties may be. As a result, 
instead of concentrating on a single indicator, we shall introduce a two parameter 
family of possible covariant indicators, which we refer to as density contrast indicators, 
Sik, hi the form 



Sir 





h ab dp 


dp 


L 


p 1 dx a 


dx b 



K 



dV, IeR, ffel\{0}. (4) 



An important feature of this family is that it may be treated as local or global 
depending upon the size of £. It also includes as special cases the indicator given 
by Tavakol & Ellis f2S|| , for which 1 = 2 and K = 1/2, and the pointwise indicator 
previously given by Bonnor [|| 

p 2 dx a dx b ' 1 ' 

In the cosmological context we might expect it to be more appropriate to normalise 
these indicators with the comoving volume. We therefore define the corresponding 
density contrast indicators per unit volume, Sikv, by 

S IK v = ?f (6) 

where V = J s dV is the comoving volume. We shall also consider dimensionless 
analogues of these indicators in section 7. 

Indicators (Q) and (JsJ) are of potential interest for a number of reasons, including 
their operational definability and hence their potential relevance to observations re- 
garding the evolution of structure in the Universe and their possible connection to the 
question of gravitational entropy 

For completeness we shall also consider the spatially covariant indicator introduced 
by Szafron & Wainwright Eq| 



l .luab. d P d P 



SW = --\ h ab -^-^r (7) 

p V dx a dx b v ' 

as well as the non-covariant indicators given for LT models by Bonnor Q 

B2= X -% (8) 
p or 

and Silk @ 

SL= r -% (9) 
p or 

where r in these expressions is the r coordinate of the LT models introduced in the 
next section. We note that some of these latter indicators have been used as measures 
of homogeneity in the past P, E6|, a question we shall return to section 9. 



■5 



In the following, in analogy with the notion of cosmological arrow which points in 
the direction of the dynamical evolution of the Universe, we shall employ the notion of 
density contrast arrow which is in the direction of the evolution of the density contrast 
indicator employed. 

The aim here is to test these families of indicators in the context of LT and Szekeres 
models in order to determine the subset of these indicators (and models) for which 
the asymptotic evolution and the evolution near the origin is monotonically increasing 
with time, i.e there is a unique density contrast arrow which points in the direction of 
increasing time. 

3. Lemaitre Tolman models in brief 

The Lemaitre-Tolman models [jl9|, ^9|, [l| (see also |l8j ) are given by 

,2 

ds 2 = -dt 2 + j^-jdr 2 + R 2 (d6 2 + sin 2 9dcf> 2 ) (10) 

where r, 6, <fi are the comoving coordinates, R — R(r, t) and / = f(r) are arbitrary C 2 
real functions such that / > — 1 and R(r,t) > 0. In this section a dot and a prime 
denote d/dt and d/dr respectively. The evolution of these models is then given by 

^ 2 = f+/ (11) 

where F = F(r) is another C 2 arbitrary real function, assumed to be positive in order 
to ensure the positivity of the gravitational mass. Equation (|ll]) can be solved for 
different values of / in the following parametric forms: 

For / < 0: 

(V ~ sin V ) = ^tlll {t - a ) (12) 

where < rj < 2n and a is a third arbitrary real function of r. 
For / > 0: 

F 

R = ^(cosh?? - 1), 

(wahri-n) = -t r (t-a) (13) 



where rj > 0. 
For / = 0: 



Ii=(^-) U-nP. (1.4) 



() 



The solutions corresponding to / > 0, / = and / < are referred to as hyperbolic, 
parabolic and elliptic, respectively. In the elliptic case, there is a recollapse to a second 
singularity, while the other two classes of models are ever-expanding. In all three cases 
the matter density can be written as 



p(r,t) = 



F 



(15) 



8ttR'R 2 ' 

Now the fact that p must be non-divergent (except on initial and final singularities) 
and positive everywhere imposes restrictions on the arbitrary functions |l6[ , with the 
positivity of the density implying that R' and F' have the same sign. 



4. Evolution of the density contrast in Lemaitre Tolman models 

For these models the indicators (||) can be written as 



Sik = 47T 



1 + / 1 



R a p 1 



K 



R 2 \R'\ 
VT+7' 



dr 



(16) 



with the time derivative 

SlK=^ J (1 + 

+ R 2 \R'\ l - 2K 



9 (R 2 lR m-2^ ] 




dp 
dr 



(17) 



In the following we consider different classes of LT models given by different types of 
/. Clearly in general p in such models depends on the functions / and F and as a 
result the integrals arising in Sik and Sikv cannot be performed in general. We shall 
therefore look at the asymptotic behaviour of these indicators and in some special 
cases we study the behaviour of the indicators for all times. 

We note that in some studies concerning the question of gravitational entropy the 
function a has been taken to be a constant or zero (see for example Bonnor j3j), in 
order to avoid the presence of white holes. Here, in line with these studies, we shall 
also take a to be a constant in all sections below, apart from section (9.2) where the 
shortcomings associated with taking a = a(r) will not affect our results. 

Parabolic LT models 

Models of this type with a = constant, reduce to the Einstein-de Sitter model [|| for 
which the density depends only on t giving Sik = = Sikv for all time. 



4-2. Hyperbolic LT models 
For large r\ we have 

R^fHt-a) 

which gives 

F' 

P 



47r/'/*(t-a) 3 
-3F' 

47r/'/i(f-a) 4 



(18) 

(19) 
(20) 
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where for positivity of p, F' must have the same sign as /'. This then gives 

4vr / ai (t-a) 3IK - SK+3 dr (21) 



J ai (t-a) 3IK - 8K dr 

Sikv ~ 2 — ! (22) 

Jlf'lfHt-afdr 

where ol\ > is purely a function of r. 

The dominant asymptotic temporal behaviour of the density contrast indicators 
described in section 2 are calculated and summarised in the column 2 of Table ([!]) and 
the conditions for Sir an d Sikv to be monotonically increasing can then be readily 
calculated and are summarised in Table (0). 



4-. 3. Elliptic LT models 

In the limit rj — > 2n, R, p and p can be written as 

2 

R, 



Mi 



3F' 



P 



27F'(—f)i ( F' 3f'\h + ^-t 



2 / 



(23) 
(24) 

(25) 



where <f> 2 (r, £) = 12 



(-/)* 



(t — a) — 7T 



&(r,t) = -12 



( EL-zL 

\F 2 / 



and 



0i (r) 



4(-/)' 



Now satisfies 2 < and 02 > and </> 2 and F' must have opposite 



signs to ensure the positivity of p. The indicators are then given by 
SjKKiK J a 2 \<f> 2 \ IK -^ K+1 dr 

Sjk w 4tt /" (jJf - yJif + l) a 2 2 |<fe| /K ^ 



+ CK2'|<^>2 



7/f- J#A'+1 



dr 



(26) 



(27) 



where a2(r,i) = l|<^ 2 |0i 



'3^4 ( 
\2 



l + f)F' 2 ( 2<f, 3 1 <p' 2 \ I 



per unit volume can be calculated to be 

_ 3 J a 2 \h\ IK -^ K+1 dr 

biKv -2 m^\d- r • 



A" 



> 0. Similarly the density contrast 



(28) 



Now with the choice of ( ^ — | 4- 



~>IK 



~>IKV- 



Whc 



3^ 
2 J) 



0, the model becomes homogeneous with 

< and our 



7^ 0, we have > and 



dt 



results are summarised in Tables ([!]) and (||). We note that the asymptotic behaviour 
of Sikv cannot be deduced in general and therefore the corresponding entries in these 
tables were derived using pointwise versions of the indicators. 



To summarise, the results of this section indicate that for both ever-expanding 
(/ > 0) and recollapsing (/ < 0) LT models, I and K can always be chosen such that 
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Indicators / > / < 



Sjk (t-a)^<-^ 4 K ~^ K+1 



Sikv (t — a) 



3IK-8K 



"2 



-2 



Bl (t — a) < m 



B2 const. cf> 2 1 



SL const. 
SW const. 



-l 



Table 1 . Asymptotic evolution of density contrast indicators given in section (0) , for 
hyperbolic and elliptic LT models. The constants in the second column are different 
for each r. 



Models Sik Sikv 



f>0 


I > 


8 
3 


1 

K 


I > 


8 
3 


f<0 


I < 


10 

3 


1 

K 


I < 


10 

3 



Table 2. Constraints on I and K in order to ensure Sik > 0; Sikv > 
asymptotically in hyperbolic and elliptic LT models. 



Sik and Sikv both grow asymptotically. However, there are special cases of interest, 
such as / = 2, for which no such intervals can be found. 

4-4- Special LT models 

So far we have studied the behaviour of the Sik and Sikv asymptotically. To partly 
extend these results to all times, we shall in this section consider some concrete ex- 
amples of LT models which have been considered in the literature. 

Parabolic examples: 

Models with a — (see e.g. those in Q and are homogeneous with trivial 

behaviour for the indicators. 

Hyperbolic examples: 

We considered a number of examples of this type given by Gibbs |]l0), Humphreys 
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et al. |TjJ, and Ribeiro (^3), the details of which are summarised in Table (||). We 
found that for all these models there exist ranges of I and K such that indicators Sik 
increase monotonically for all time. In particular, we found that the condition for all 
time monotonicity with K = 1/2 and K = 1 are given by / € [1, +oo[ and / £ [2, +oo[ 
respectively. 



References 


F(r) 




i 


Humphreys et al. 


F = y 


f = 


r 2 


Humphreys et al. 


F = ±r 3 


f = 


r 3 


Gibbs 


F = F tanh r 3 


f = 


/o sinh 2 r 


Ribeiro 


F = F rP 


f = 


sinh 2 r 



Table 3. Examples of hyperbolic LT models, where Fo,p and /o are positive 
constants. 



Elliptic examples: 

We considered examples of this type, given by Bonnor Q and Hellaby & Lake |Q, 
details of which are summarised in Table (|^). Again we found that for all these models 
there exist ranges of / and K such that indicators Sik are monotonic for all time. In 
particular, we found that the condition for all time monotonicity with K = 1/2 and 
K = 1 are given by / S ]0, 1] and I G ]0,2] respectively. For values of K = 1/2 and 
/ outside the range / < 4/3, e.g. I = 2, we find that Sji increases in the expanding 
phase while decreasing after a certain time (which depends on r) in the contracting 
phase and tending to zero as the second singularity approaches. 



References f(r) 
Bonnor F = F r 3 f = -/ 0T ^r 

Hellaby & Lake F = F ^ f = -/o^r 

Table 4. Examples of elliptic LT models, where _Fo and fo ^ 1 are positive real 
constants and m, n are integers such that m> n. 



To summarise, we have found that for all these concrete examples there exists 
values of / and K such that indicators Sik increase monotonically for all time. 
However, the allowed intervals of / and K are in general narrowed relative to those 
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obtained by the asymptotic considerations. Finally the indicators Sikv m &U these 
models lead to non elementary integrals. 

5. Szekeres models in brief 

The Szekeres metric is given by [£7| [n] 

ds 2 = -dt 2 + R 2 e 2u {dx 2 + dy 2 ) + R 2 H 2 W 2 dz 2 (29) 

where W — W(z) and v — v(x, y, z) are functions to be specified within each Szekeres 
class and R = R(z,t) obeys the evolution equation 

M 

R 2 = -k + 2—, (30) 

where M = M(z) is a positive arbitrary function for the class I models and a positive 
constant for class II models, defined below. The function H is given by 

H = A- p + f + -p-f- (31) 

where functions A — A(x,y,z), H, R and W are assumed to be positive; p + and 
/3_ are functions of z and /+ and /_ are functions of z and t, corresponding to the 
growing and decaying modes of the solutions X of the equation 

•■ 2RX 3MX , x 

X + — -^T = °- ^ 
The density for these models is given by 

3MA , s 

p{ - x,y,z,t ' = ^Rm- (33) 

The solutions to (|30| ) are given by 

dhirj) 
R = M — j-^- 
dr\ 

t - T{z) = Mh{rf) (34) 



where 



rj — sin 77 (k = +1) < 77 < 2tt 
Hv) — { smhr/-i] (k = -1) < i] (35) 
\rf (k = 0) < 77, 

which in turn allows the solutions to ( |32| ) to be written as 

( QMR- 1 ^- |?7cotij?)-l k = +l 
f + (r])=i 6MR-\1 - ^cothiyy) + 1 k = -1 (36) 

I ToV 2 ' ' k = 

( 6MR- 1 cot ±r) k = l 
f_(r))=) 6MJ?- x cothir7 k = -1 (37) 

[ 24?T 3 fc = 0. 

The Szekeres models are divided into two classes, depending upon whether d ^ e z ^ 7^ 
or 9 ^ ) = 0. The functions in the metric take different forms for each class and for 
completeness are summarised in the Appendix. 
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6. Evolution of the density contrast in Szekeres models 



In these models the density contrast indicators are given by 



SiK — 



1 



dp 
dx 



1 



e^R 2 p I \dy 



H 2 W 2 R 2 p I \d 



dp 



2v 1 
K 



dp 



R 3 e 2u HWdxdydz 



(38) 



We note that the singularities in these models are given not only by R = but also 
by H = 0, which define the so called shell crossing singularities. For example, choos- 
ing /3-f > eventually results in a shell crossing singularity at a finite time given by 
A = /?+/+ (see jllj for a detailed discussion). Here in order to avoid shell cross- 
ing singularities, we either assume (3+ < in all cases, or alternatively f3+ > and 
f3+(z) < A(x,y,z), in the k — — 1 case. 

We recall that for class II models, T = constant. In the following we shall, in line 
with the work of Bonnor [H , also make this assumption in the case of class I models 
(implying /3_ = 0) in order to make the initial singularity simultaneous and hence 
avoid white holes. 



We consider the two Szekeres classes in turn and in each case we study the 
three subclasses referred to as hyperbolic, parabolic and elliptic, corresponding to 
the Gaussian curvatures k = —1,0,-1-1 respectively. 



6.1. Evolution of the density contrast in class I models 

Assuming /3 + = in this class makes p = p(t)f, which implies Sik = = Sikv f° r 
all time. We shall therefore assume (3+ to be non-zero and possibly z dependent. The 
contribution of /3_ would be important near the initial and final singularities but is 
irrelevant for the asymptotic behaviour of parabolic and hyperbolic models. 



Parabolic class I Szekeres models: 



For the models of this type with T = constant, the density depends only on t 
which trivially gives Sik = = Sikv for all time. 



Hyperbolic class I Szekeres models: 



For large r\ we have 



resulting in 



R »i-T (39) 
+ l (40) 



3MA 

PW 4tt{A -f3+)(t -T) 3 [ ' 



f In fact, the Szekeres models, with T = constant, reduce to FLRW iff /?+ = /3_ = 



12 

. _ -9 MA 

P ~ 4ar(A-p+)(t-T)* ^ ' 

which have the same t dependence as the corresponding LT models. Using (38) we 
obtain 

Sir ~ J a 4 (t - T) 3KI - SK+3 dxdydz (43) 

^ Ja4{t-T) 3KI - 8K+3 dxdydz 
IKV ~ J e 2 »W(A - /?+)(£ - Tf dxdydz ( ' 

where is a positive function of x, y, z. 

The dominant asymptotic temporal behaviour of the density contrast indicators 
for these models and the conditions for their monotonic behaviour with time are sum- 
marised in Tables (||) and (j^). 

Elliptic class I Szekeres models: 

Using (|33"|), together with the following approximations 



R « (45) 



U « ^ (46) 



where — (^jf - — 27r), gives 
6A 

(47) 

— 6A 

(48) 

which has again the same temporal dependence as for the corresponding closed LT 
models. Using (38) we obtain 

Sik ~ J ct5^pi K 3 K+1 dxdydz (49) 

2 J a 5 ip 2 dxdydz 
IKV ~ 9 / Me^Wp+^dxdydz ' [ ' 

where 



a 6 (x,y,z,t) = \MWe 2v -K^ 



K 

96 2 48~ J 6~5 A-'e- 2 " \ ( ( dA\ 2 j_ ( dA\2 j_ e 2l >/ 

A/ir-' 2 - r (7r/3 + ) 2 --f 



('Mr _L (dA\2 _r_ e ^3 



and as 77 -> 2?r, i -> 2-kM + T and -02 -> -2ttM'. 



Our results are summarised in Tables (|5|) and Q. We note that the asymptotic 
behaviour of Sikv cannot be deduced in general for these models and therefore the 
corresponding entries in these tables were derived using pointwise versions of the 
indicators. 
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6.2. Evolution of the density contrast in class II models 

Recall that in this class T and M are constants in all cases. Again if (3 + — /3_ = we 
recover the homogeneous models giving Sik = = Sikv- In what follows we shall 
therefore consider the general cases with (3+ — f3+(z) and /3_ = (3-{z). 

Parabolic class II Szekeres models: 

The asymptotic evolution of p in this case is given by 

P* — r (51) 

M-P+)(t-T)i 

P~ -■ (52) 

97T/?+(t-T)T 

Using (38) we obtain 

Sik » J a e (t~T)i KI -^ K +ldxdydz (53) 

s ^ 2 ja 6 (t-T)i KI -^ K +idxdydz 
IKV ~ 9 Je^\(3+\MW(t~T)idxdydz 

where a§ is a positive function of x, y and z only. Our results are shown in Tables ([s]) 
and (|). 

Hyperbolic class II Szekeres models: 

Here we use the approximations for R, H, /+ and /_ already shown for class I. The 
asymptotic evolution of p and dp/dt are in this case given by 

P ~ «{A-fi + ){t-T)s (55) 
. -18MA 

'"^-/MCt-D* (56) 

Using (38) we obtain 

Sik « / «r(t - T) 3IK - SK+3 dxdydz (57) 

fa 7 (t-Tr IK -* K + 3 dxdydz 
IKV ~ / e 2 -lU(^l - /3+)(t - Tfdxdydz 1 j 

where 07 is a positive function of x, y and z. Our results are depicted in Tables 
and (|). 

Elliptic class II Szekeres models: 

The asymptotic evolution of /? and its time derivative are in this case is given by 

9 ~ ^M 2 (7r/3 + -/3_)^2 { ' 



-6A 



P-n7^ ^72 (60) 
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where ip^(z, t) = (rjj- — 27r). Here we impose the restriction 7r/3 + 
the positivity of p and by (38) we have 



S 



IK 



KI- 



SlKV 



as I i>2 
/ "sIV^ 



'-K+1 



dxdydz 



/3_ < to ensure 



(61) 



KI-MKA 



dxdydz 



27 / e 2v {ir(3 + - P_)Mip 2 dxdydz 



(62) 



where as is a positive function independent of t. The asymptotic behaviours of the 
indicators for these models are identical to the corresponding models in class I (except 
in the parabolic case) and the conditions for their monotonic evolution with time are 
summarised in Tables (||) and (^J). Again the asymptotic behaviour of Sikv cannot be 
deduced in general for these models and therefore the corresponding entries in these 
Tables were derived using pointwise versions of the indicators. 





Class II 


Classes 


I & II 




Indicators 


k = 


k = -1 


k = 


-- +i 


Sir 


(t 


- T)i IK ~^ K + 2 


(t — T) 3IK ~ 8K+3 


IK- 

v 2 




Sikv 


(t 


-T)l IK ~f K 


(t — j i y iIK - 8K 


IK- 


fK 


Bl 


(t 


-T)-i 


(t-T)- a 






SW 


(t 


-T)-i 


const. 


4 





Table 5. Asymptotic evolution of a number of density contrast indicators for the 
class I and class II Szekercs models. Columns 3 and 4 represent the behaviours for the 
hyperbolic and elliptic models which are identical for both classes, while the second 
column represents the behaviour for class II models. 



To summarise, the results of this section indicate that for both ever-expanding 
(k = and k — —1) and recollapsing (k = +1) Szekeres models, we can always choose 
I and K such that the Sm and Sjkv & re asymptotically increasing both separately 
and simultaneously. However, for some cases of interest, such as e.g. I = 2, no such 
interval can be found for which both sets of indicators simultaneously have this prop- 
erty. Also as can be seen from Table (jj|) different indicators can, for different values 
of I and K, give different predictions concerning the asymptotic homogenisation of 
these models. 

We also note that the similarity between the results of this section (for k = +1 
and k = — 1) and those for LT models (for / > and / < 0) is partially due to 
the fact that the density functions p has the same time dependence in both of these 
sub-families of LT and Szekeres models. A nice way of seeing this, as was pointed out 
to us by van Elst, is that for the above dust models the evolution equations for the 
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Sik 


SlKV 


k = 






J ^ 2 ft 


1 ^ 2 


k = -1 


J ^ 3 ft 


i> s 


J ^ 3 ft' 


J ^ 3 


k = +1 


J ^ 3 if 


J < io 

^3 


J ^ 3 ft" 


/< io 

^3 



Table 6. Constraints on 7 and K in order to ensure Sik > 0; Sikv > 
asymptotically in class I and class II Szekeres models. 



density, expansion, shear and electric Weyl curvature constitute a closed dynamical 
system which is identical for both models (see 0). This, however, does not necessarily 
imply that our indicators should also have identical time evolutions for all times for 
both models, since they also include h ab and dV in their definitions. It turns out that 
asymptotically they are the same in the cases considered here. 



7. Evolution of dimensionless indicators 

As they stand, indicators (||) are not dimensionless in general. To make them so, we 
shall also briefly consider their dimensionless analogues given by 

Sikl = (63) 
where L is a real number which depends on I and K thus 

L = —IK — 2K + 1. (64) 
3 

Asymptotic behaviour of Sikl for LT and Szekeres models are summarised in Table 
(0). This demonstrates that there are still intervals for I and K (I € ]2,4[, K € R\{0}) 
such that these dimensionless indicators asymptotically increase with time. In this way 
the results of the previous sections remain qualitatively unchanged. 



8. Behaviour near the initial singularity 

So far we have studied the asymptotic behaviour of these models at late times. 
To further understand the possible monotonic behaviour of these indicators, it is 
of interest to study their behaviour near the initial singularities. Our results are 
summarised in Table (j^) and the constraints on I and K in order to ensure the non- 
divergence of Sik, Sikv and Sikl near the singularities are given in Table (|9|). 

As can be seen from Table (JsJ) , apart from the case of Szekeres II with decaying 
modes (J3— ^ 0), indicators have the same behaviour as we approach the origin for 
all other models. On the other hand, in presence of decaying modes, the constraints 
on I and K necessary to ensure initial non-divergence and final monotonic increase 
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Models 



Tolman 



Szekeres I 



Szckcres II 



Parabolic (t-T)% IK -i K 
Hyperbolic (t - a) IK - 2K (t-T) IK ' 2K {t-T) H< - 2K 



Elliptic 



klK-^K 



IP, 



klK-^K 



1>: 



klK-±K 



Table 7. Asymptotic behaviour of Sxkl f° r LT and Szekeres models. In the elliptic 
cases pointwise versions of the indicators were used. 



are disjoint, ic, there are no intervals of I and K such that these conditions are 
simultaneously satisfied. This in turn seems to imply that there is an incompatibility 
between the presence of decaying modes (see also [| [TT| ) and a unique density 
contrast arrow. 



Indicators Tolman 



Szekeres I 



Szekeres II 

W- + o) 



Szekeres II 

W- = o) 



SlK 

SlKV 

SlKL 



(t - a ) 2IK - iK + 2 (t - T ) 2IK - iK + 2 (t - T) 

\2IK-iK (+ rp\IK-^K 



(t - a) 2IK ^ K 
(t-a)i IK 



{t-Tf 
(t-T)i IK 



(t-T)? IK -i K 



(t — j i ^ 2IK -^ K + 2 
(t - T) 2IK - AK 
{t-T)i IK 



Table 8. The behaviour of the indicators Sik, Sikv an d Sikl in LT and Szekeres 
models near the initial singularity. 



9. Some consequences 

In the following we shall briefly discuss some of the conceptual issues that the above 
considerations give rise to. 

9.1. Expanding versus recollapsing phases 

An interesting question regarding cosmological evolution is whether the expanding 
and recollapsing evolutionary phases do (or should) possess the same characteristics 
in terms of density contrast indicators and in particular whether or not the rise of 
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Lemaitre-Tolman 
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Szekeres I 
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> 2- 
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I > 2 


/if > 


Szekeres II = 0) 
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> 2- 


1 

K 


I > 2 


/if > 


Szekeres II ^ 0) 


I 


> IS _ 

^ 3 


1 

if 


^3 


/ > 4 



Table 9. Constraints on i and A" in order to ensure the non-divergence and 
monotonic increase of Sik, Sikv ar >d Sjkl near the origin of time. 



structuration is expected to increase monotonically throughout the history of the 
Universe independently of the cosmological direction of its evolution, in other words, 
whether the density contrast and the cosmological arrows have the same sign. 

The answer to this question is of potential importance not only with regard to the 
issue of structure formation but also in connection with debates concerning the nature 
of gravitational entropy and its likely behaviour in expanding and contracting phases 
of the Universe (see e.g. Hawking Q, Page plj and also references in H). We 
note that the underlying intuitions in such debates seem to have come mainly from 
considerations of the Friedman-Lemaitre-Robertson- Walker (FLRW) models together 
with some studies involving small inhomogeneous perturbations about such models 
p5| . Our studies here, though classical, can therefore provide potentially useful com- 
plementary intuitions by considering exact inhomogeneous models. In particular, our 
results in the previous sections raise a number of points that can potentially inform 
this debate, amongst them: 

(i) Indicators in recollapsing models: We have found that different indicators behave 
differently in expanding and recollapsing phases of cosmological evolution, with 
some that remain monotonically increasing in both phases and others which 
change their signs in the expanding and recollapsing phases. In particular we 
find that there is a larger class of indicators that remain monotonic in the ever- 
expanding models and in this way we could say that it is more difficult to maintain 
monotonic growth in recollapsing phases. 

(ii) Spatial dependence of turning time in recollapsing models: As opposed to the 
closed FLRW models, the turning time in inhomogeneous models can be position 
dependent. For example, in the LT models with a = 0, the turning time is given 
by t — 7rF/(2/2), which in general depends on r. In this sense there are epochs 
over which the Universe possesses a multiplicity of cosmological arrows. 

This raises the interesting question of whether there can be observationally viable 
inhomogeneous cosmological models which allow ranges of epochs and neighbour- 
hoods, such that over these epochs different neighbourhoods (labelled Eat in the 
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definition of the indicators (4), (6) and (p3|)) can give rise to different local density 
contrast arrows within the same model. 

We note that in the ever-expanding models this problem is less likely to arise, 
which raises the question of whether this can be taken as an argument in favour of 
ever-expanding models which more easily allow uni-directionality in their density 
contrast arrows. 

9.2. Connection to homogeneity 

Another question is what is the connection between homogeneity and the behaviour of 
density contrast indicators? We start by recalling a result of Spero & Szafron |^{| ac- 
cording to which Szekeres models are spatially homogeneous (with Sik = = Sikv) 
if the density in these models is a function of t only. Therefore in these models 
Sik = ==>■ homogeneity. Now in general Sik = may not imply homogene- 
ity, which raises the interesting question of what is the set of inhomogeneous models 
which satisfy this property? A related question has been considered by van Elst 
who points out some of the difficulties involved. We hope to return to this question 
in future. 

Another point in this regard is that different indicators (even covariant ones) can 
make contradictory statements about whether or not a model homogenises asymptot- 
ically. Here to illustrate this point we look at the following examples. 

First, we consider the parabolic LT model studied by Bonnor Q in the form 

ds 2 = -dt 2 + (t- a) 4 / 3 ( (l + 3(f^y) df 2 + f 2 dcA (65) 

where a = a(r). Employing the indicator B2, Bonnor deduced that for fixed f this 
model approaches homogeneity, as t — > oo, irrespective of its initial conditions. This 
indicator is not covariant. But even for covariant indicators (such as Sik) one can 
find ranges of I and K (namely I > ^ — and its complement) which give opposite 
conclusions regarding asymptotic homogenisation. 

As another example, we consider the hyperbolic LT models studied by Bonnor Q 
and given by F = bf 3 / 2 : with b = constant. According to the Bl, B2, SL and SW 
indicators these models approach homogeneity, whereas there are ranges of values of 
I and K (namely I > ^ — - 1 ?) for which our measures Sik increase monotonically. 

We also take the hyperbolic class of LT models with a = 0, which were studied in 
the subsection (4.2) above. As can be seen from Table (|J), there are ranges of / and 



K for which the asymptotic behaviour of both Sik and Sikv can be the opposite to 
the Bl indicator. 

Finally, It is also of interest to compare our results with the studies of asymptotic 
behaviour for Szekeres models. An interesting result in this connection is due to Goode 
& Wainwright [jll), according to which Szekeres models with k = 0,-1 and [3 + = 
become homogeneous asymptotically, a result that was also obtained for the class II 
models by Bonnor & Tomimura [o . We recall that for the models with simultaneous 
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initial singularity considered here, the assumption of f3 + = reduces the class I models 
to that of FLRW. For the class II models (with /3_ ^ 0), the asymptotic behaviour of 
our indicators depends on the choice of the indices I and K. In particular, to ensure 
the asymptotic increase of the indicators Sik , Sikv, Sikl for k — and k = — 1, it 
is necessary and sufficient to have I > Q — j?, I > I > 5 and 7>^ — i,7>^, 
I > 4 respectively. 

10. Conclusions 

An important assumption in gravitational physics concerns the tendency of 
gravitational systems to become increasingly lumpy with time. Here, we have tried to 
study this possibility in the context of inhomogeneous models of LT and Szekeres, by 
using a number of two parameter families of density contrast indicators. 



3.0 - 




Figure 1. The region between the curves gives the ranges of values of / and K 
such that Sik is increasing both near the singularity and asymptotically for all the 
models considered here. The squares show the special examples of indicators with 
I = 1,K = 1/2 and 1 = 2,K = 1. 

Given the arbitrary functions in these models, we have only been able to 
establish conditions for monotonicity of our indicators as we approach the origin and 
asymptotically. Even though these are necessary but not sufficient conditions for 
monotonicity for all times, nevertheless our studies illustrate some general points that 
seem to be supported by our all time study of a number of special models. Our results 
show: 

(i) Different density contrast indicators can behave differently even for the same 
model. We find there is a larger class of indicators that grow monotonically 
for ever-expanding models than for the recollapsing ones. In particular, in the 
absence of decaying modes = 0), we find that indicators exist which grow 
monotonically with time for all the models considered here. Figure (1) gives a 
brief summary of our results by depicting the ranges of I and K such that Sik 
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grow monotonically near the origin as well as asymptotically for all the models 
considered here. An example of a special indicator that lies in this range is given 
by K = 1, 1 = 2 (a non pointwise version of Bl). However, the indicator given by 
K = 1/2,1 = 2 (which is linear in the derivatives of density (28)) does not satisfy 
this property. 

(ii) If decaying modes exist (i.e. /3_ ^ 0), we find no such indicators which grow 
monotonically with time for all the models considered here. Recalling a theorem 
of Goode and Wainwright [[n] , namely that (3- = is the necessary and sufficient 
condition for the initial data (singularity) to be FLRW-like, our results seem 
to imply that the presence of monotonicity in the evolution of density contrast 
indicators considered here is directly related to the nature of initial data. This 
is of potential relevance in the recent debates regarding the nature of arrow of 
time (to the extent that such an arrow can be identified with the density contrast 
arrow) in ever expanding and recollapsing models (T 

(iii) Our considerations seem to indicate that the notion of asymptotic homogenisation 
as deduced from density contrast indicators may not be unique. 

(iv) The possible spatial dependence of turning points in inhomogeneous models can 
lead to multiplicity of local density contrast arrows at certain epochs, which could 
have consequences regarding the corresponding cosmological arrow(s) of time. 

Finally we note that given the potential operational definability of our indicators, 
it is interesting that the overall approach to homogeneity may not necessarily be 
accompanied by a decrease in the density contrast, as measured by the different 
indicators. In this way different covariant density contrast indicators may give different 
insights as to what may be observed asymptotically in such inhomogeneous models. 
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11. Appendix 

The different forms of the functions in Szekeres class I and II models 



Class I 

R = R(z,t) 
/± = /±(M) 

T = T(z) (66) 
M = M(z) 

e v = f[a{x 2 + y 2 ) + 2bx + 2cy + d}- 1 



2i, m. 



n 
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where functions in the metric are subjected to the conditions: 



ad _ b 2 _ c 2 = 1 £ 

T^(z) 2 = (e-fc/ 2 )- 1 
/?+(*) = -kfMzQM)- 1 
— fT z (GM)~ 1 



e = 0,±1 



(67) 



The functions a, b, c, d and / are all functions of z that are only required to satisfy 
equations (67). 



Class II 



R = R(t) 
f± = f±(t) 

T = const. (68) 



M = const 
e v = [1 
W = 1 



e» = [l + ±k(x 2 + y 2 )}- 




in this case a, b, c, /?+ and (3- are arbitrary functions of z. The curvature is given by 
k. 
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